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Abstract: We consider the zero heat conductivity limit to a contact discontinuity for the 
mono-dimensional full compressible Navier-Stokes-Fourier system. The method is based on the 
relative entropy method, and do not assume any smallness conditions on the discontinuity, nor 
on the BV norm of the initial data. It is proved that for any viscosity > 0, the solution of the 
compressible Navier-Stokes-Fourier system (with well prepared initial value) converges, when 
the heat conductivity k tends to zero, to the contact discontinuity solution to the corresponding 
Euler system. We obtain the decay rate . It implies that the heat conductivity dominates the 
dissipation in the regime of the limit to a contact discontinuity. This is the first result, based 
on the relative entropy, of an asymptotic limit to a discontinuous solutions for a system. 


1 Introduction 

We study the zero dissipation limit of the solution to the Navier-Stokes-Fourier system of a 
compressible heat-conducting gas in the Lagrangian coordinates in one-dimensional space M: 

Tt-Vx = 0, 

^ Vt+piT,^)x = 

/ / nx / / ax N / X /X 

(e(T, 'll) -b —)t + (p(r, i^)v)x = k(—) x + - )x 

'-2 T T 

where the functions T{x,t) > > 0 represent the specihc volume, velocity and 

the absolute temperature of the gas, respectively, p = p{T,'d), e = e(r,-d) is the pressure and 
the internal energy satisfying the second law of thermodynamics, > 0 is the viscosity constant 
and K > 0 is the coefficient of heat conduction. 

*The work of A. F. Vasseur was partially supported by the NSF Grant DMS 1209420. E-mail: 
vasseur@math.utexas.edu. 
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If K = = 0 in (|l.ll) . then the corresponding Euler system reads as 


Tt - rza; = 0, 

vt+p{T,'d)x = 0, 

(e(T,'(9) + y)^ + iPv)x = 0. 

For the Euler equations (11.21) . it is well-known that there are three basic wave patterns, i.e., 
shock and rarefaction waves in the genuinely nonlinear helds and contact discontinuities in the 
linearly degenerate field. If we consider the Riemann problem of compressible Euler system (II.2h 
with the following Riemann initial data 




x < 0, 

(tr, vr, ^r) X > 0, 


(1.3) 


then the above Riemann problem admits a contact discontinuity solution 


{T,v,'d){x,t) = I 

provided that 

VL = VR := V, 

The relative entropy method was first introduced by Dafermos m and Diperna m to prove 
the stability and uniqueness of Lipschitzian solutions to the conservation laws. The method 
was used by Leger in [28] to show the stability of shocks in (up to a shift) in the scalar case 
(see also [T] for an extension to L^, 1 < p < oo). The method has been extended to the system 
case in [29| for extreme shocks. The method is purely nonlinear, and allow to work without 
smallness assumptions. The relative entropy method is also an effective method for the study 
of asymptotic limits. One of the first use of the method in this context is due to Yau [38] for 
the hydrodynamic limit of Ginzburg-Landau models. Since then, there have been many works 
in this context, see ima nia laEoi [32] etc. and the survey paper [36], although they are 
all concerning the limit to a smooth (Lipschitz) limit function. Recently, the relative entropy 
method has been successfully applied in to prove the vanishing viscosity limit of the viscous 
scalar conservation laws to shocks m, and the contraction of viscous shock profiles to the 
scalar viscous conservation laws [22]. The theory of stability of discontinuous solutions, based 
on the relative entropy has been reformulated in [331 [23] in terms of contraction, up to a shift. 

This paper is the first attempt to use this theory for the study of asymptotic limits to a 
discontinuous solution, in the case of systems. In the present paper, we are concerned with 
the vanishing physical viscosity limit of compressible Navier-Stokes-Fourier system to the 
above contact discontinuity solution (|I.4p to the Euler system (|1.2p . with arbitrarily large wave 
strength, using the relative entropy method. The contact discontinuities of the Euler system 
enjoys a contraction property of a special inhomogeneous pseudo-distance, which does not need 
a shift (see m)- This property is one of the key of our proof. 


iTL,VL,'&L) 

(rR^VR^'&R) 


X < vt, 
X > ut. 


(1.4) 


r)r = nr> •= n 
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For the zero dissipation limit of conservation laws, there are many results using BV-theory, 
see Bianchini-Bressan[6] etc. and the limit to basic wave patterns to the compressible Navier- 
Stokes equations and Boltzmann equation by the elementary energy method, one can refer to 
[Ml EU [H Eol ED and the references therein. 

In the present paper, we consider the state equation as in Feireisl m 

=Pe{T) +Pth{T,'&) (1.5) 

with Pe{T) being the elastic pressure. The thermal pressure pth{T,'&) takes the forms as 

Pth{r,'&) = t9p^(t). 

Both functions Pe('r) and p^(t) satisfy 

Pe(T), Pair) > 0, Pe(r), p'^ij) < 0. (1.6) 


Correspondingly, the internal energy takes the form as 

e(r,??) = Pe(r) + Q(i?), (1.7) 

where Pe{T) is the elastic potential dehned by 

-Pe(T) = -y Peiy)dy, 

and (5('!?) = f'^Cy(z)dz is the thermal energy satisfying 

— Cy{d), 

where Cy{'d) is a positive and smooth function of d. By the second law of thermodynamic 


dds = de + pdr, 


we can calculate the entropy s as 

s{T,d) = J ^'’^^^ dz - (1.8) 

with ^ 

P^{t) = - j P^{T)dT. 

Set U = (t,v, e(T,d} + being the solution to the compressible Navier-Stokes equations 
(11.11) and U = (r, v, e(f, d) + ^)* being the contact discontinuity (II.4p to the compressible Euler 
equations (|1.2p . As usual, the associated relative entropy is defined by 

s(PlU) := s{U) - s{U) - ds{U) • {U - U), 
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with s being the entropy in (11.81) . Correspondingly, the relative entropy-flux is 

q{U;U) ■.= -ds{U)-{f{U)-fiU)), 

with the flux f{U) = {—v,p{T,'&),p{T,'d)vY. 

Dehne the entropy functional as 


Sit) = - I ?? siU\U)dx. 


(1.9) 


Direct calculations show that for the state equations in (II.5j) and (ini), one has the explicit 
formula for (|1.9I) 

S{t) = — J -d siU\U)dx = J [ — ^ 'S(^) + ^ s(C7) -|- ?? ds\jj^fy ■ [U — U)] dx 


[ ^^dz -'dP^f) + ^ J ^^dz 


Ppir — t) — v{v — v) + (E — E) 


dx 


( 1 . 10 ) 


-!?(P^(r) - P^(r) -hp^(r)(T - r)) {Peir) - Peir) -FPe('r)(r - r)) 


+ {Qi^) - Qi^) - ^ [_ 


^C,iz).. 1 


dz) + -iv - 

■& z A 


dx. 


Now we state the main result as follows. 


Theorem 1.1. For k> 0, let (t,.^, "d^) a weak solution to dnD such that siUf^\U) lies in 

L^(0,T;L^(R)) where = (t^," i^K) and U is the contact discontinuity to the Euler system 
m) and s is defined in dUD. Assume further that the pressure p satisfies that 


Pe[r) r^T 

with 7 > 2, 

clS T - y 0“t“; 

(1.11) 

P^{t) ~ 

with 7 > 2, 

cis T — y 

(1.12) 


Then if the initial perturbation satisfies that 

Sio) = O(V^), 

and 

s{Uo^\Uo) G L“(D) (1.13) 

with D C R being any neighborhood of x = 0. Then for any ia > 0, there exists kq > 0 such that 
if 0 < K < kq, then it holds that 

S{t) = 0(-v/k), Vt G [0,T] with any fixed T > 0. 

Remark 1.2. There is no need of uniform-in-n bound on the weak solution 11^, even though 
its existence is still open now. However, if both the initial volume and initial temperature have 
the lower and upper bound, then existence of global classical solution can be proved by following 
the similar procedure as in Kazhikhov-Shelukhin Jl^ . 
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Remark 1.3. By (|1.10l) . it can he seen that \\vk — However, for — t and 

one only can get the convergence in entropy norm as in (fTTol) and can not derive the 
convergence in — norm unless one has the uniform lower and upper bound for both and 
a priorly. 

Remark 1.4. Theorem [LJ\ shows that the solution of compressible Navier-Stokes system can 
converge to contact discontinuity solution to the corresponding Euler system with the decay 
rate uniformly for any viscosity v > t), which implies something surprising that the heat 
conductivity k dominate the dissipation in the regime of limit of contact discontinuity as k tends 
to zero. Note that the wave strength of contact discontinuity can be arbitrarily large and the 
initial perturbation can contain the possible vacuum states formally. 

Remark 1.5. For the state equation, the pressure p can contain a large class of gases, such as 

. \ Rd 

P{T,d) = - hPe{T) 

T 

with R > 0 being the gas constant and the elastic pressure Pe{T) satisfying the condition (ll.lip . 


2 Proof of Theorem 1.1 


In this section, we will prove Theorem 1.1 by relative entropy method. We will first define the 
cut-off entropy functional by suitably choosing the cut-off function r/(x) based on the underlying 
structure of contact discontinuity. The cut-off function r]{x) is chosen to satisfy 


r]{x) 


1, X < —1, 

smooth and decreasing, — 1 < x < 1, 

0, X > 1, 


and 

rj'{—x) = 7?^(x). 

For definiteness, we can choose t/(x) as follows. 


7?(x) 


'1, X < -1, 

1 o 3 1 

< 4 ^ “ 4 ^+ 2 ’ 

0, X > 1. 


Obviously, we have 


h'{x) 



3 

4’ 


-1 < X < 1, 

otherwise, 


satisfying r]' < 0, supp rj' = [—1,1] and the symmetry in (12.2p . 
Define the cut-off entropy functional as 


X 

£i{t) = - r,{-)dL s{U\UL)dx 

J —oo ^ 



r]{-^)dR s{U\UR)dx 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


(2.4) 


(2.5) 
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where e > 0 represents the layer-width for the contact discontinuity in the vanishing dissipation 
limit and is to be determined as e = ^/K in the sequel. The functional £ i { t ) and the cut-off 
function r](x) both play an important role in our analysis. 

In fact, the following proof is independent of the state equations in (jl.5p and (jl.7p until the 
estimation (|2.14|) . that is, all the proofs before (I2.14p are valid to the generic gas satisfying the 
second law of thermodynamic. By the second law of thermodynamics, it holds that 


'dds 


de + pdr = d{E 


—) +pdT 


dE — vdv + pdr, 


where the total energy E 


e(T,i9) + 


thus the entropy s satishes the equation 


,'dT. dl v: 




( 2 . 6 ) 


By direct calculations, it holds that 
d £ i { t ) 


dt 




[ ( h (-) + h (--)) (pVx - — - Px ) + — — - { Pv ) x ) dx . 

./—/-V-, ^ S S ' 7~ 7~ 7” / 


(2.7) 


Note that the third term on the right hand side of (j2.7|) vanishes due to the integration by part 
and the symmetry of the p'{x). In fact, 


1-00 ^ - Px ) + - { pv ) x)dx 

/ rY> /y \ / 


X. , X 


V - V ) x - - ( P - P ) x ) 

+ “ {Pv-pv)a^dx 


X. , X 


I +00 

-OO 


J —OO o \ 6 c / \ 


_ Vx , . Vx VVx , _ 

V — v) — VV - {p — p) + K - h U - {pv — pv) ] |lc 

r T T 

Vx dx VVx _N 1 , 

V — V) — vv - \P — P) + Vi -h V - [pv — pv) dx 

T T T 


( 2 . 8 ) 


= 0 . 


By integrating (j2.7p with respect to t over [0, t] with t £ [0, T], wow we arrive at 

+ /-OO, ^ ’'<-7>'’«) ("tI + 

= ^ 1 ( 0 )-^ j (vi^)dL + p{-^)dR^K{^)jxdt. 


(2.9) 
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Note that the terms on the right hand side of (11.101) are non-negative under the condition (11.61) . 
Then it holds that 


Si{t) = - 


f-l-OO 




r >+00 


'n{-)'&L s{U\UL)dx - 


' e 


s{U\UR)dx 


rO 


,x. 


f-l-oo 


X, 


> / J](-)[-dL s(U'll7L)]dx + / 7]( -)[-dR s([/lU'R)]dx 

J —oo ^ ./ 0 ^ 

1 1 r + oo 

[-'dLs{U\UL)]dx + - [-'dR s{U\UR)]dx 

^ J-oo ^ Jo 

= ^ s([/iu}]dx = 


( 2 . 10 ) 


Substituting the above inequalities (|1.10p - (l2.10|) into (12.91) . one has 


1 r+co 

- / i9(P^(t) - P^(f) +P^(f}(T - f)) + (Feir) - Pe{f) + Peif){T - f)) 

^ J —OO ^ ^ 

+ [Q{d) - Q{d) - d !_ dz) + ]r{v- v)"^ {x,t)dx 

% 


J-d z 

1 p r + °° ,;2 

L 1-0. 


<Ti(0) - 


( 2 . 11 ) 


Now we treat the terms on the right hand side of (j2.1ip . First, for the initial values satisfying 
(jl.lSp . it holds that 


/ -|-00 /* ~|“00 

r){-)dL smUL)dx - / ii{--)dR s{Uo\UR)dx 

^ -OO ^ J— oo ^ 

= [ 'ni-)[-'dL s{Uo\UL)]dx + [ ri{-)[-'dL s{UQ\UL)]dx 

J — OO ^ Jo ^ 

X X 

+ v{ - )[-dRs{Uo\UR)]dx+ r]{ - )[-dR s{Uo\UR)]dx 

Jo £ J-£ £ 

P + OO f£ 


< 


-do s(UolUo)]dx + [ r]{-)[-dL s{Uo\UL)]dx 
Jo £ 


+ / vi—)[-'dR s{Uo\UR)]dx 


< 


J —e ^ 0 

T(0)-|-C' f r]{—)dx-\-C [ r]{ - )dx<£{0) + C£ 

Jo £ J-£ £ 


( 2 . 12 ) 


for some generic positive constant C provided that e = y/K is suitably small. Now it remains 
to control the last term on the right hand side of (12.111) . In fact, by integration by part and 
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Cauchy inequality, it holds that 


//■ 


rt r+oo / T \ 

+ vi - )^R 

Jo J-oo V e 

rt r+oo J X, ^ \ KtJx , , 



0 J —CO 


< ^min{7?L,'i?ij} 


i:l 


■t re 

0 J-e 



v'{^){^L-^R)^^dxdt 

^lo J_^Wi^)f^dxdt 


K iJn 


£ rd 

%4«+ (t)L-tlnr » 


min{t?L,'i?ij} e 


< 


1 /■* r d^ 

-mm{dL,dR} J J K^p^dxdt + C{dL,dR) 


(2.13) 


K 



1 


dxdt 


0 J —£ 


where the positive constant C{dL,dR) = ig • With the notation 1{...} being the char¬ 
acteristic function of the set {• • • } C R X [0, T] and Ci being suitably large positive constant to 
be determined, it holds that 


J ^dxdt = J ^{l^i^Ci} + '^{^>Ci})dxdt 


-l^i^u^^dxdt. 


(2.14) 


Now the state equations in (jl.5p and (jl.7p as in [13] are crucially used to control the second 
term on the right hand of (|2.14p . Moreover, if the pressure function Peix) satisfies (II.lip or 
P^{t) satisfies pi.l2p . then 

Peir) or P^{t) ~ 7 > 2, as r ^ 0 -|- . 


Therefore, we can choose the suitably large positive constant Ci snch that 

1 

^ ^{P-ffix) - P^X) + p^{f){T - f)) + {Peir) - Pe{f) + Peif){T - f))^ {^X^i} < 
therefore, one has 


pt poo 

; c / / [d{p^{x) - p^{f) + p^{f){T - f)) 

Jo J —oo 


+ {Peix) - Pe{x) +Pe{x){x - x))]dxdt. 


Substituting p2.16p into p2.14p and then into (I2.13P yields that 

1 /■* P d'^ K P r°° 

<-mm{dL,dR}J J K^dxdt + C{dL,dR)-^ J J [??(P^(r)- 

+P^{f){x - r)) {Peix) - Pe{f) +Peix){ 


-PAf) 

T — f))~\dxdt. 


(2.16) 


(2.17) 

















The combination of (12.111) . (I2.12P and (12.171) implies that 

ft /■+00 „,2 


11 /■* 7 +°° 

-£{t) + + 

t/ 0 — CXD 


< £ (0) Ce C I £-\- 


Set 


F{t) = f £{t)dt, 
Jo 

then applying Gronwall inequality to (|2.18l) implies that 

F(t) <C2(e + 


T(0) + c-^t 


for some positive constant C 2 . Furthermore, one has 


£(t) < C* T(0) 


T(0) + £\ c-^t 


e 


(2.18) 


(2.19) 


( 2 . 20 ) 


Choosing the contact discontinuity layer e = ^/K and the initial perturbation T(0) = O^^/k) as 
in Theorem 1.1, one has 

£{t) = 0{^), Vts[0,r], Vz/>0, 

which completed the proof of Theorem 1.1. 
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